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I. INTRODUCTION 



Chiral perturbation theory (ChPT) [l|, 0, 0] is the effective field theory of the strong 
interactions at low energies (for an introduction see, e.g., j3, H, @])- It relies on a perturbative 
expansion in terms of small parameters q/A, where q denotes a quantity like the pion mass 
or external momenta that are small relative to the scale A, which for ChPT is expected to be 
of the size 1 GeV. One of the essential ingredients of ChPT is a consistent power counting, 
which assigns a chiral order D to each Feynman diagram for the process in question and 
which predicts that diagrams of higher orders are suppressed. Assuming the coefficients of 
the perturbative expansion to be of natural size one would expect contributions at order D+l 
to be suppressed by a factor q/A compared to contributions at order D. For q of the order of 
the pion mass and A m 1 GeV, this corresponds to a correction of about 20%. In the mesonic 
sector of ChPT this rough estimate seems accurate, however, the situation is less clear for 
the baryonic sector. While for example the chiral expansion of the nucleon mass shows a 
good convergence behavior, the nucleon axial coupling receives large contributions from 
higher-order terms M]. Further examples include the electromagnetic form factors of the 
nucleon (see, e.g., [1, 0]), which only describe the data for very low values of momentum 
transfer. For some of these quantities higher-order contributions clearly play an important 
role. The description of the nucleon form factors can be improved by the inclusion of vector 
mesons as explicit degrees of freedom, which corresponds to the resummation of higher-order 
contributions [1, ^ . 

The convergence properties of baryon chiral perturbation theory (BChPT) are also of 
great importance for lattice QCD. Due to numerical costs, present lattice calculations still 
require pion masses larger than the physical one, and results obtained on the lattice have 
to be extrapolated to the physical point. ChPT as an expansion in the pion mass is the 
appropriate tool to perform such extrapolations, which again poses the question for which 
values of small parameters the ChPT expansion gives reliable predictions. 

There are several renormalization schemes for manifestly Lorentz-invariant BChPT at 
the one-loop level that result in a consistent power counting while preserving all relevant 
symmetries [l(J 11, 12^13, 14, 15]. The most commonly used of these is the infrared (IR) 



regularization of Ref. [llj. All these renormalization schemes have in common that there is a 
relation between the chiral order and the loop expansion, so that the investigation of higher 
chiral orders requires the evaluation of multi-loop diagrams. In Ref. [l6j a reformulated 
version of the IR regularization has been introduced that is also applicable to multi-loop 
diagrams [l7[. Reference [l8j contains a different generalization of IR regularization to 
two-loop diagrams. 

In this paper we describe the renormalization procedure for two-loop integrals in mani- 
festly Lorentz-invariant BChPT within the reformulated version of infrared regularization of 
Ref. 0. The results of a calculation of the nucleon mass up to and including order 0(q 6 ) 



have been reported in Ref. [19]. To the best of our knowledge, this is the first complete two- 
loop BChPT calculation in a manifestly Lorentz-invariant framework. Here, we describe the 
details of the calculation. In particular we show that the renormalization procedure pre- 
serves all relevant symmetries and that renormalized two-loop diagrams obey the standard 
power counting rules. A calculation of the nucleon mass to order Q(q 5 ) was performed in 
Ref. [20] in the framework of HBChPT (see, e.g., jU, H3]), and Ref. [23| contains the leading 
non-analytic contributions to the axial-vector coupling g& at two-loop order obtained from 
renormalization group techniques. 
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This paper is organized as follows. In Section [TT] we review the main features of infrared 
renormalization at the one- loop level that are essential for the following. Section ITO1 contains 
a brief overview over the general aspects of the renormalization of two-loop diagrams. The 
infrared renormalization of products of one-loop integrals is described in Sec. IIV} while the 
discussion of genuine two-loop integrals follows in Sec. EJ An application of these methods 
can be found in Sec. Wl\ which is followed by a summary. Explicit expressions for the 
appearing integrals can be found in the appendix. 



II. INFRARED REGULARIZATION OF ONE-LOOP INTEGRALS 



The method of infrared regularization was developed as a manifestly Lorentz 



invariant renormalization scheme preserving all relevant symmetries. It results in diagrams 
obeying the standard power counting (see Sec. I VI p . Infrared regularization is based on di- 
mensional regularization and the analytic properties of loop integrals, and in its original 
formulation is applicable to one-loop integrals containing pion and nucleon propagators in 
the one-nucleon sector of ChPT. In Ref. 16] a different formulation of infrared regulariza- 



tion was presented which reproduces the results of the original formulation up to arbitrary 
order. The advantage of the new formulation is that it can also be applied to multi-loop 
diagrams and diagrams containing additional degrees of freedom 0, [lg, [13, 0, 25]. Since 



IR regularization can, in fact, be viewed as a renormalization scheme we also refer to it 
as infrared renormalization. We briefly describe those features of the renormalization of 
one-loop integrals which are important for the renormalization of two-loop integrals. 
Denote a general one-loop integral containing pion and nucleon propagators by 

/d n k 1 
(27r) n a x ---a m b x ---bi 

where a, = (k + qi) 2 — M 2 + i0 + and bj = (k + pj) 2 — m 2 +i0 + are related to pion and nucleon 
propagators, respectively, and n = 4 + 2e is the number of space-time dimensions. Infrared 
renormalization consists of splitting the integral into an infrared singular part I n ...N--- and 
an infrared regular part R-^.-.n- -, 

II ::■ X- = I :: X- + R. X- , (2) 

or for short 

H = I + R. (3) 

The advantage of splitting the original integral into two parts is that the infrared singular 
part Itt -n - satisfies the power counting, while R n ...N -- contains terms that violate the power 
counting. In addition, the infrared singular and infrared regular parts differ in their analytic 
properties. For noninteger n the expansion of in small quantities results in only 

noninteger powers of these variables, while R-^...^... only contains analytic contributions. In 
the formulation of Ref. [16| one obtains the infrared regular part R^...^... by reducing H V ... N ... 



to an integral over Schwinger or Feynman parameters, expanding the resulting expression in 
small quantities such as pion masses or small momenta, and interchanging summation and 
integration. 

As an example consider the integral 

*M0, -P)=1 J ^y ^ 2 _ M2+ iQ+ ] p _ p )2 - m 2 + i0 +] ■ (4) 
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To apply the reformulated version of IR renormalization we combine the two propagators 
using 

1 f 1 dz 
ab 



,„ \{l-z)a + zbY 1 
and perform the integration over the loop momentum k, resulting in 

HMO, -P) = r(2 - n/2) f dz [C(z)] in ^- 2 



(47T)n/2 



(5) 







where C(z) 



2^2 



m z 



(p 2 - m 2 )(l - z)z + M 2 (l - z) - i0 + . Next, we expand [C(z)] 



(n/2) -2 



in p — m and M and interchange summation and integration. This generates the chiral 
expansion of the infrared regular part R, which is given by 



R 



m n " 4 r(2-n/2) 
' (47r)"/ 2 (n-3) 



p 2 — m? (n — 6) {p 2 — m 2 ) 2 (n — 3)M 2 



2m 2 



4m 4 (n - 5) 



2m 2 (n — 5) 



(6) 



and which coincides with the expansion of R given in Ref. [111 ]. 

Symmetries introduce relations among various Green functions of the theory, called Ward- 
Fradkin-Takahashi identities (Ward identities for short) j26l . 27 . 28| . Expressions contain- 
ing the integrals H satisfy the Ward identities, 1 since they are derived from an invariant 
Lagrangian and dimensional regularization does not violate the symmetries. Since J, for 
noninteger n, only contains nonanalytic terms, while R consists of analytic contributions 
only, each part has to satisfy the Ward identities separately in order for the sum H = I + R 
not to violate any symmetry. 

Both the infrared regular and the infrared singular parts contain additional divergences 
not present in the original integral H. Since these additional divergences do not appear in 
H, they have to cancel in the sum of / + R = H. This means that 



R 



add 



radd 



(7) 



The e expansion of H is given by 

H uv 



H 



+ # (0) + 0(e) 



e 

J-fUV jadd j^add 

+ + + I + R, 



(8) 



where H uv /e denotes the ultraviolet divergence of H, H^ ' refers to the terms independent of 
e, and we have explicitly shown the additional divergences in the second line. In BChPT the 
renormalization can be performed in a two-step process. First, all divergences are absorbed, 
and then additional finite terms are subtracted. In the standard approach the divergences 
are absorbed using the MS scheme. In this scheme one subtracts the quantity 



327T 2 



- ln(47r) + 7s - 1 



In the following we use the phrase that integrals satisfy the Ward identities, by which we mean that 
expressions containing these integrals satisfy the Ward identities. 



4 



where = — r'(l), and sets the appearing t'Hooft parameter [29| fi = m, where m is 
the nucleon mass in the chiral limit. Here, in order to simplify the calculation, we apply 
minimal subtraction (MS) with a t'Hooft parameter fi, absorbing only terms proportional to 

e , and then set fi = 7^72 e^ - . This is completely equivalent to the standard approach. 
The infrared renormalized expression H r of the integral H is defined as its finite infrared 
singular term, 

W = J, (9) 

which satisfies the power counting since all terms violating it are contained in R. One of the 
fundamental properties used in the construction of the effective Lagrangian is the invariance 
under symmetries of the underlying theory. It is therefore of utmost importance that these 
symmetries are not violated at any step in the calculations. We now show that the definition 
of the renormalized integral H r of Eq. (Q satisfies this requirement [ll| . The original integral 
H is obtained from a chirally symmetric Lagrangian using dimensional regularization, which 
preserves all symmetries. Therefore expressions containing H satisfy the Ward identities; 
and in particular their e expansions satisfy the Ward identities order by order. As explained 
above, R satisfies the Ward identities separately from I. This also means that the Ward 
identities are satisfied order by order in the e expansion of R and J, respectively. Therefore 
the identification of the renormalized integral H r as H r = I does not violate any symmetry 
constraints. Since the sum of additional divergences cancels, the term which is subtracted 
from H is given by 

H uv - 

R=—— + R. (10) 
With Eq. ([7]) and the definition of Eq. (jTUl) we can write 

H = I + R. (11) 
Within the framework of dimensional regularization, the dimensional counting analysis 



of Ref. [30J provides a method to obtain expansions of loop integrals in small parameters. 
This method is described in detail in Appendix |A] Here we show how the infrared regular 
and infrared singular parts of the integral H are related to the different terms obtained from 
this method. Using dimensional counting, H is written as 2 

H = G 1 + G 2 . (12) 

For G\ we simply expand the integrand in M and interchange summation and integration. G2 
is obtained by rescaling the integration variable kv-^^k and then expanding the integrand 
with subsequent interchange of summation and integration. For p 2 = m? the method of 
obtaining G\ is the same as the one used to determine the expansion of the infrared-regular 
part R. It follows that 

G\ = ^ R n = R, (13) 

n 

while G2 gives the chiral expansion of the infrared singular term /, 

G 2 = J2 J n, (14) 



2 Note that the notation used here differs from the one found in the Appendix to avoid confusion with terms 
in the e expansion of H . 
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(a) 



(b) 



(c) 



(d) 



FIG. 1: Two-loop diagram with corresponding subdiagram and counterterm diagram. 



where R n and I n are the terms in the chiral expansion of the infrared regular and infrared 
singular parts, respectively. It should be noted that the expansion of I does not always 



converge in the entire low-energy region [11|. For the integrals considered in the calculation 
of the nucleon mass, however, the expansion of I converges. The identification of G\ and 
G*2 with the infrared regular and infrared singular parts, respectively, is used below to show 
that the renormaiization process in the two-loop sector does not violate the considered 
symmetries. 



III. GENERAL FEATURES OF THE RENORMAIIZATION OF TWO-LOOP IN- 
TEGRALS 

We give a brief description of the general renormaiization procedure for two-loop integrals 
before presenting details of the IR renormaiization. The discussion follows Ref. [311 ]. 

At the two-loop level integrals not only contain overall UV divergences, but can also 
contain sub divergences for the case where one integration momentum is fixed while the 
other one goes to infinity. As an example consider the two- loop diagram of Fig. [TJ (a). It 
contains one- loop subdiagrams, shown in Fig. [TJ (b). The renormaiization of subdiagrams 
requires vertices as shown in Fig. [Tj (c), which are of order h. 3 At order h 2 these vertices 
appear in so-called counterterm diagrams as the one shown in Fig. [Tj (d). When the sum 
of the original diagram and the one-loop counterterm diagrams, Fig. [Tj (a) and twice the 
contribution from Fig. [TJ (d), respectively, is considered, the remaining divergence is local 
and can be absorbed by counterterms. In order to renormalize a two-loop diagram one has 
to take into account all corresponding one-loop counterterm diagrams. 

In our calculation we encounter two general types of two-loop integrals. The first type 
can be directly written as the product of two one-loop integrals, while this decomposition is 
not possible for the second type. 

IV. INFRARED RENORMALIZATION OF PRODUCTS OF ONE-LOOP INTE- 
GRALS 

Consider the product of two one-loop integrals, 

H = H a H b . (15) 



3 Here, the power of h denotes the order in the loop expansion. 
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if is a two-loop integral and the result of a dimensional counting analysis reads (see App. lA 2|) 



H = H^ + + H^ + 

= F 1 + F 2 + F 3 + F 4 , (16) 

where F±, F 2 + F 3 , and F 4 satisfy the Ward identities separately due to different analytic 
structures, i.e. different overall powers of M in n dimensions. Using Eq. ([3]), H can also be 
expressed as 

H = I ah + hRb + R a h + R a Rb, (17) 

where again I a I b , I a Rb + Rah, and R a Rb satisfy the Ward identities individually. 

To renormalize the integral H we need to add the contributions of (renormalized) coun- 
terterm integrals. The vertex used in the counterterm integral is determined by standard IR 
renormalization of a one-loop subintegral. In a one-loop calculation we do not have to con- 
sider terms proportional to e for the subtraction terms since, at the end of the calculation, 
the limit e —>■ is taken. At the two-loop level, however, the subtraction terms are multi- 
plied with terms proportional to e _1 from the second loop integration. Therefore the choice 
whether or not to include the terms proportional to e in one-loop subtraction terms results 
in different finite contributions in the two-loop integrals. In addition to the UV divergences 
and the terms proportional to e° we choose the subtraction terms for one-loop integrals to 
contain all positive powers of e, 4 

ttUV 

R = ^— + R^ + eR^ + --- . (18) 

This choice is crucial for the preservation of the relevant symmetries as is discussed in the 
following. H contains two subintegrals, H a and Hf,. The expressions for the unrenormalized 
counterterm integrals then read 

— R a H b — RbH a , (19) 

which themselves need to be renormalized applying IR renormalization. The Hi (i = a, b) 
are one-loop integrals from which we would subtract the term Ri in a one-loop calculation, 
excluding the additional divergences. However, the term Rj multiplying Hi contains terms 
with positive powers of e, so that in the product of Rj and Ri we get finite terms from the 
additional divergences in Ri (see Eq. ([7j)). These would not be removed if we chose the 
subtraction term to be RjRi. Instead we define the subtraction term for the product RjHi 
to be 

JjUV jjadd JDadd 

— RjRi + ' y +#Ti-, (20) 

i.e. we subtract all finite terms stemming from the additional divergences in Ri but do not 
subtract the additional divergences themselves. This is analogous to the one-loop sector, 
where we do not subtract the additional divergences in the infrared regular part either (see 
Eqs. (ED and (HOD). 

We now show that this renormalization procedure for the counterterm integrals does not 
violate the Ward identities. We know that the subtraction terms S for one-loop integrals do 
not violate the Ward identities and result in a modification of the coupling constants and 



4 In a calculation at two-loop order it is sufficient to include terms proportional to e. Higher powers of e 
are required for the generalization to multi-loop diagrams. 
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fields in the Lagrangian. The counterterm integrals are then calculated with the help of this 
new Lagrangian which means that the term 

-SHi (21) 

also respects all Ward identities. Hi is a one-loop integral and Eq. (f2"Tj) can be written as 

-Sh-SRi (22) 

where —Sli and —SRi satisfy the Ward identities separately. In particular, the Ward 
identities are satisfied term by term in an expansion in e for S Jj and S Ri, respectively. The 
expansion for S h is given by 



' odiv \ / jadd \ qdiv jadd i 

g j / |_ gf in | / Zl |_ j/* n | i _|_ _ div j fin _|_ gfinjadi 1 



e I \ e " / e 2 e 



+ ••■• (23) 



Suppose we choose the finite part of the counterterm to vanish, 5 

gfin = q 

In this case we can see that the term proportional to e -1 is given by 

Igdivjfin (24) 

It has to satisfy the Ward identities since for this choice of S it is the only term proportional 
to e -1 in the e expansion of S Ii. By changing the renormalization scheme to also include 
finite terms in the subtraction terms, the product in Eq. (|24p does not change, but we 
obtain the more general expression of Eq. f[2"3"j) . Considering the term proportional to e _1 
and keeping in mind that Eq. f l2~21) respects the Ward identities we now see that 

- S fin If d (25) 

satisfies the Ward identities separately. Since the additional divergences have to cancel in 
the sum of Ii and Ri it follows that If dd = —Ri dd and 

- ~ S fin Rf d (26) 

does not violate any symmetry constraints. Using the fact that S Ri respects all symmetries 
and choosing the subtraction term S to be Rj (which only contains UV divergences), 

C p ndiv £jUV qfin 

O — j"i — 3 ' — 3 ' 

it follows that 

rrUV -nadd -padd 

Rf^l_ (27) 



5 In baryonic ChPT this would result in terms violating the power counting. So far we are only concerned 
with the symmetries of the theory, which are conserved for S' m = 0. The issue of power counting is 
addressed below. 
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satisfies the Ward identities and therefore also our prescription for the subtraction terms of 
the counterterm diagrams of Eq. (|2Q|) satisfies the Ward identities. 

Using the above method the sum of the original expression and the renormalized coun- 
terterm integrals gives 

JjUV TDadd TDadd TDadd ZjUV 

U U p rr I p p a p(0) p jt , p p n b 



e 2 e 



TDadd 



e 



JjUV JD^add _ JDadd 

I ah + -^a-Rfe + -^fe-Ra + -Ra-Rfe — -Ra-^fi 2 — ^ ^1 — " Rbh 

e z e 

TDadd TjUV TDadd 
K " U b 6(0) K 

— /t^ 



e 2 " b e 



_ JjUV JDadd _|_ jDadd jjUV _ jDadd 

h + l a (Rb — Rb) + h(R a — R a ) - — - — n — - — R a — 

e 2 e 

R b — - h R a Rb ■ (28) 



The difference between Ri and Ri is only given by the additional divergences R® dd / e, resulting 
in 

jadd \ / j add \ / jadd \ TDadd / jadd \ TDadd 

I + V) ( h + V) + + V) *r + ( h + V) 

ZjUV TDadd i TDadd ttJJV TDadd TDadd 

_ H a R b + R a H b _ R{ , } R^_ _ R( o)R^_ + RaRb (29) 
Using If dd = —R? dd we obtain 

jadd jadd JjUV TDadd i TDadd jtUV TDadd JDadd 

Iah - - Q Rb b ~ Ra° ] ^- ~ RT— + RaRb + O(c). (30) 

Expanding R a Rb in e and simplifying the resulting expression gives 

jadd jadd JDadd JDadd ttUV ttUV ttUV p(°) _i_ p(°) ttUV 

U b -±4- + ^ LJ T ! -+ 2 + — & +(RaRb) {0) + O(e), (31) 

e z e 

where (R a Rb)^ stands for the terms proportional to e° in the product R a Rb- Using again 
jadd _ _Radd we gee a ]] t erms containing the additional divergences vanish, 

ttUV ttUV ttUV 6(0) , 6(0) ttUV 

Iah + + a b e a " + (RaRb) {0) + 0(e) . (32) 

The term R a Rb satisfies the Ward identities, in particular each term in the e expansion of 
R a Rb does so individually. This means that we can subtract the finite part of R a Rb by 
a counterterm. The terms proportional to e~ 2 and e _1 stem from the UV divergences in 
H a and Hb- These terms also satisfy the Ward identities individually and are absorbed in 
counterterms. As desired, the renormalized result for the product of two one-loop integrals 
including the counterterm integrals is then simply the product of the renormalized one-loop 
integrals, 

(H a H b ) r = Iah- (33) 
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Besides respecting all symmetries the renormalization prescription must also result in 
a proper power counting for renormalized integrals. The chiral order of a product of two 
integrals is the sum of the individual orders. For a one-loop integral the infrared singular 
part I satisfies the power counting. Therefore the result of Eq. (1531) also satisfies power 
counting. 



V. INFRARED RENORMALIZATION OF TWO-LOOP INTEGRALS RELE- 
VANT TO THE NUCLEON MASS CALCULATION 

In this section we describe the renormalization procedure for two-loop integrals that do 
not directly factorize into the product of two one-loop integrals. We follow the general 
method presented in Ref. 17j], but give more details. First we show how the proper renor- 



malization of two-loop integrals and of the corresponding counterterm integrals preserves 
the underlying symmetries. Next we describe a simplified formalism to arrive at the same 
results while greatly reducing the calculational difficulties. 

A. General method 

Denote a general two- loop integral contributing to the nucleon mass by H 2 , 

f f d n k\d n k 2 1 

Ho(a, b, c, d, e\n) = // — — — : , . — , (34) 

J J (2n) 2n A a B b C c D d E e ' v ; 

where 

A = k\ — M 2 + i0 + , 
B = kl-M 2 + iO + , 
C = k 2 1 +2p-k 1 +i0 + , 
D = k\ + 2p- k 2 + i0 + , 

E = k 2 l +2p-k 1 + 2k 1 -k 2 + 2p-k 2 + kl + i0 + . (35) 

Using a dimensional counting analysis we can write H 2 (a, b, c, d,e\n) as 6 

H 2 = F 1 + F 2 + F 3 + F 4 . (36) 

Fi is obtained by simply expanding the integrand in M and interchanging summation and 
integration. For F 2 we rescale the first loop momentum k\ by 

k\ i— > — ki, (37) 
m 

expand the resulting integrand in M, and interchange summation and integration. is 
obtained analogously to F 2 , only that instead of k\ the second loop momentum k 2 is rescaled, 

M 

k 2 ^— k 2 . (38) 
m 



6 For brevity we employ the notation H2 for #2(0, b, c, d,e\n) in the following discussion. 
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Finally F4 is denned as the result from simultaneously rescaling both loop momenta, 

ki h-> — fci, fc 2 ^ — ^2, (39) 
m m 

and expanding the integrand with subsequent interchange of summation and integration. 
F±, F 2 + F 3 , and F4 separately satisfy the Ward identities due to different overall factors of 
M. This is analogous to the one-loop sector, where the infrared singular and infrared regular 
parts separately satisfy the Ward identities, since the infrared singular part is nonanalytic in 
small quantities for noninteger n, while the infrared regular term is analytic. As in the one- 
loop case the interchange of summation and integration generates additional divergences not 
present in H2 in each of the terms Fi, F2 + F3, and F4. Again, these additional divergences 
cancel in the sum of all terms. 

In addition to the two-loop integral we also need to determine the corresponding subin- 
tegrals. To identify the first sub integral we consider the k\ integration in H2, 

f dnkl 1 

subi ~ J (2<r A-c-E- ■ { ' 

This is a one-loop integral which is renormalized using "standard" infrared renormalization. 
The infrared regular part -R S u&i of this integral is obtained by expanding the integrand in M 
and interchanging summation and integration. The only term in Eq. (HOI depending on M 
is A. Symbolically we write 

_ f dnkl 1 

Rsu bl - 2^ J , (41) 

where underlined expressions are understood as an expansion in M. R SU bi contains additional 
divergences, and we define R su b 1 as R su bi without these divergences, 7 

j^add 

Rsuh = RsvJbi • (42) 

As in the definition of Eq. (jSJ), R su b 1 again contains all terms of positive powers of e. Since 
H subl is a standard one-loop integral, R su b 1 will satisfy the Ward identities and can be 
absorbed in counterterms of the Lagrangian. 

Using these counterterms as a vertex we obtain a counterterm integral of the form 

HcTi = ~ S Wr Rsubl ' (43) 

Hctx is generated by a Lagrangian that is consistent with the considered symmetries. There- 
fore, H C t± satisfies the Ward identities. Inserting Eqs. (|4T!) and (|42l) we rewrite FL C t x as 

r d n k 2 v r d n h 1 r d n k 2 Rf u d bl 1 

CTl J (2n) n ^J (2n) n A a B b C c D d E e J (2tt)« e B b D d ' 1 ' 

Equation fT44|) still needs to be renormalized. After the ki integration has been performed, 
Eq. (1441) is a one-loop integral and standard infrared renormalization can be used. To 



7 Note that for the integrals of interest here, the UV divergence is included in the infrared regular part R. 
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obtain the infrared singular part Ict x we rescale k 2 i— > — k 2 , expand in M, and interchange 
summation and integration. Symbolically we write 



CTl ^7 (2tt)" ^7 (27r)™A^B 6 C^D d £ e 



V / f - l E add 1 (4^ 

J ( 27r )« 5^ ' 



where double-underlined quantities are first rescaled and then expanded. Note that R^L 
can also depend on k 2 through the denominator E in Eq. (1401) . Since Ict\ is obtained 
from a one-loop integral that satisfies the Ward identities through the standard infrared 
renormalization process, it will itself satisfy the Ward identities. The infrared renormalized 
expression for the counterterm integral, 

jadd 

I CTl = Ict x - ^f- , (46) 

then also satisfies the Ward identities. Note that 1^ itself contains terms proportional to 
e _1 , since it stems from the one-loop counterterm for the subintegral, but we choose not 
to include any terms proportional to positive powers of e. This means that e" 1 /^ only 

contains terms proportional to e~ 2 and e~ l . The expression for Ia\ therefore does not 
contain any divergent terms stemming from additional divergences. 

We now show how Ict y is related to the term F 3 of Eq. (|36|) . As explained above, F 3 is 
obtained by rescaling k 2 , expanding the resulting integrand and interchanging summation 
and integration. In the above notation this would correspond to 

F y 1 (47) 

^JJ {2ix) 2n A a g?C c D^E^' V ; 

Comparing with the first term in Eq. (|45j) we see that the integrands in both cases are 
expanded in the same way. Therefore, when adding the counterterm diagram Iqt x to H 2 it 
cancels parts of F 3 . The difference between Iqti and F 3 is that in F% the terms stemming from 
the additional divergences R^f (including finite terms) as well as the additional divergences 
IcT./e that are proportional to e -2 and e -1 are not subtracted. As pointed out above, the 
original integral H2 only contains UV divergences, therefore the additional divergences cancel 
in the sum Fi + F 2 + F 3 + F 4 . Apart from these contributions, the terms remaining in the 
sum Iqtx + ^3 are the finite contributions stemming from the additional divergences in 
Rsubx ■ Since in F3 the variable k 2 is rescaled before expanding while the k\ variable remains 
unchanged, F3 can be considered as a sum of products of infrared singular and infrared 
regular terms, which we symbolically write as 

F 3 = J2 R ^- ( 48 ) 

In this notation the remaining finite terms are ^2,R1 dd I^\ where e~ l R1 dd is the additional 
divergence of R\ and l!p is the part of I 2 proportional to e. 

The second subdiagram can be calculated analogously, and is related to the term F 2 in 
Eq. <m. 
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Taking the above considerations into account we obtain for the sum of the original integral 
H 2 and the corresponding counterterm integrals 



H 2 + Icti + Ict 2 — F\ + F 2 + F 3 + F4 + l C T x + Ict 2 

= Fx + F 4 + R a i dd l2 ] + Yl R 2 ddl i ] 

= Fi + F± - ti* 1 ® ~ Yl J 2 dd A (1) , (49) 

where Fi indicates that the additional divergences are excluded. 

The expression in Eq. ( 1491) satisfies the Ward identities since each term in the sum on the 
left side of the first line does so individually. F\ separately satisfies the Ward identities, in 
particular this is the case for each term in its e expansion. This means that we can subtract 
the finite part of F\ by an overall counterterm without violating the symmetries. Since the 
remaining UV divergences also satisfy the Ward identities, absorbing them in an overall 
counterterm does not violate the symmetries. The result for the renormalized two-loop 
diagram is then 

HI = F 4 -J2 A ocM 4 1} - J 2^ J i (1) • (50) 

Since all subtractions preserve the symmetries if£ will satisfy the Ward identities. 

So far we have subtracted pole parts in the epsilon expansion. Following [ll| we choose 
to absorb the combination 



(47T) 



2 



1 1 

- (log(4vr) - 7b + 1) 



n-4 2 



(51) 



instead, which is achieved by simply replacing the original t'Hooft parameter fi by 



(52) 

(see also App. ID]) . 

F4 is obtained by rescaling both h\ and hi and satisfies the power counting. Since the 
terms result from the rescaling of ki, the product has the same analytic structure in 
M as F4, and therefore satisfies the power counting. This means that also the renormalized 
integral if£ = F 4 — I^ dd I^ — I^I^ obeys the power counting. 



B. Simplified method 

In the previous subsection we have established the concept of infrared renormalization of 
two-loop integrals. The procedure outlined above is quite involved when applied to actual 
calculations of physical processes. Therefore, we now describe a simpler method of obtaining 
the renormalized expression if£ which, however, is only applicable to integrals with a single 
small scale. This is the case for the calculation of the nucleon mass, whereas e.g. the nucleon 
form factors contain the momentum transfer as an additional small quantity. 

Instead of calculating the sub integrals of the original integral H 2 , consider just the terms 
in F4. F4 itself is a sum of two-loop integrals. Each two-loop integral contains one-loop 
subintegrals, which correspond to performing only one loop integration while keeping the 
other one fixed. These subintegrals contain divergences, resulting in divergent as well as 
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finite contributions when the second loop integration is performed. In addition to the subin- 
tegral contributions, contains finite parts and additional divergences originating in the 
interchange of summation and integration when generating We can symbolically write 
F 4 as 

Fa F a Fa Fa 



F 4 + ^— + + — F " 2 + ~ F t • (53) 



Here, the finite parts of F4 are denoted by F4 to distinguish them from F4 in Eq. (l4*9|) . The 
bar notation is also used for the divergent terms F% ' 2 and F% ,l to show that these are not 
the complete divergent expressions for F4, but only the additional divergences of order e~ 2 
and e _1 , respectively. The terms e -^p^ ub ^ dm denote the divergences of the subintegral with 

k • 

respect to the integration over k i} while F A 3 stands for the remaining second integration of 
the counterterm integral. Note that the divergent part of the first loop integration over ki 
in general depends on the second loop momentum kj. This dependence is included in the 
expression F* 3 . 

We now show how the different parts in Eq. (1531) are related to expressions in F 2 and F 3 
and then describe the simplified renormalization method. F4 is obtained from the original 
integral H by rescaling k\ and k 2 , expanding the resulting integrand in M and interchanging 
summation and integration. For the denominators of Eq. fl35l) the rescaling results in 

k 2 M i + t0 + ^ (MYftf m 2 + i0+), 
\m J 

k 2 -M 2 + iO + » (^j\k 2 2 -m 2 + t0 + ), 



ki + 2p ■ k x + i0 + i-> — — ki + 2p ■ k x + iO H 

\m J \m 

9 , . f M\ ( M , 9 



ki + 2p ■ k 2 + i0 + (-»• — — k 2 + 2p- k 2 + iO 

\m J \m 

k 2 + 2p-k x + 2k x -k 2 + 2p-k 2 + kl + i0 + ^ 

— — k\ + 2p ■ h + 2— h ■ k 2 + 2p ■ k 2 + — k\ + i0~ 
\m J \m m m 

After the interchange of summation and integration one can perform the substitution ki 1— > 
ki to bring the denominators k\ — m 2 + i0 + and k\ — m 2 + i0 + back into the form A 
and B, respectively. The result can be interpreted as obtained from the original integral by 
leaving A and B unchanged and expanding C in k 2 , D in k 2 , and E in k\ + 2k\ ■ k 2 + k 2 , 
respectively. Symbolically 

d n k x d n k 2 1 



F 4 ~J2 



x 



(2n) 2n [k 2 — M 2 + iO+] a [k 2 — M 2 + i0 + } b [ki 2 + 2p-kx + i0+} c 

1 

\k 2 2 + 2p-k 2 + tO+Ylh 2 + 2p-h + 2h -k 2 + 2p-k 2 + k 2 2 + iO+Y ' 



(54) 



where we have used the underlined notation to mark terms that we have expanded in. 

The divergent parts of the k\ subintegral stem from the integration region k% — > 00. 
They can be generated by further expanding each term in F4 in inverse powers of k\. This 
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corresponds to an expansion in positive powers of M for the first denominator and in positive 
powers of 2p ■ k 2 in the resulting last propagator, 

F Sub u div ^ ^ ffd n k 2 d n k t 1 



e ^JJ (2tt) 2 " [k 2 - M + i0+} a [k 2 - M 2 + i0+} b 

1 

V 

[h 2 +2p-h + iO+Ylh 2 + 2p-k 2 + i0+] d 

x \ (55) 

[h 2 + 2p-h + 2fci ■ k 2 + 2p ■ k 2 + k 2 2 + io+y ' 

We see that the expression for 2 is of the form 

pk 2 ST~^ f d n k 2 ffiu\- k^k v 2 k\ ■ ■ ■ , . 

4 ~ 2^ J (2tt)« [k 2 — M 2 + i0+} b [2p ■ k 2 + iO+Y+n ' {00) 

where f^ v \ - denotes the coefficients that result from the expansion in Eq. (I55I) . 

Next we show that F* 2 is related to terms in F$. F% is generated from the original integral 
if 2 by rescaling k 2 , expanding the resulting integrand, and interchanging summation and 
integration. After the substitution k 2 i— ► || k 2 and using the above notation we write 



+ c 



x 



-E 



d n k x d n k 2 1 
(2vr) 2ri jjfe| — M 2 + zO+] a [fc| — M 2 + iO+} b [kf + 2p-h + i0 

1 

[fc 2 2 + 2p ■ k 2 + i0+] d \k 2 + 2p-h + 2k, -k 2 + 2p-k 2 + k 2 2 + i0+] e 
d n hd n k 2 [k 2 ]^ [2fci • fc 2 p' 5 

(2tt) 2 ™ [A; 2 + i0+] a +^ [/c 2 + 2p • + i0+] c+e +^ [A; 2 — M 2 + i0+] 



i, 



1 

X [2p ■ k 2 + iQ+] d +i3 ' (57) 

We see that F% is the sum of products of one-loop (tensorial) integrals. As explained above 
these products of one-loop integrals are in fact products of infrared singular and infrared 
regular parts of integrals (see Eq. 



^3 = ^2 R^ 2 ' 

and the expressions for I 2 are given by 

j ST^ f d n k 2 k^k^k^ ■ ■ ■ , . 

2 ~ ^ J (2vr)« [k 2 — M 2 + iQ+} b [2p ■ k 2 + i0+] d +^ ' 1 ' 

Considering the k 2 integrals of Eqs. ( l55l) and ( IBTl) one sees that one has expanded in the 
same quantities. While the ordering of the expansions as well as the interchanges of summa- 
tion and integration are different, the two expansions are equivalent. Therefore, comparing 
Eqs. ( 1561) and (|58[) . one finds that for each term in F^ 2 there is a corresponding term in J 2 , 
or symbolically 

Ft 2 =l 2 . (59) 
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An analogous analysis for the second subintegral gives 



Ft 1 = h ■ (60) 

As a next step we show that the divergences of the F 4 subintegrals are related to the 
additional divergences of the integrals Ri in F 2 and F 3 . From Eq. (1551) we see that the 
divergent part of the ki subintegral is given by integrals of the type 



d n ki 



(2ir) n [kf — M 2 + i0+} a [ki 2 + 2p-h + i0+]< 
1 

\k\ + 2p ■ ki + 2fci ■ k 2 + 2p ■ k 2 + k 2 2 + iO+Y 

»l~v _ _ _ 



d n ki k^k 



(61) 



^ J (2tt)™ [kf — M 2 + i0 + ] a [2p ■ k x + i0+] c + e +^ ' 
The infrared regular integrals R\ in Eq. (|57p read 

r r^X^ [ ^ 1 ^i^i ; ; ; f 62 ) 



{2i\) n [fcf + i0+] a+mi [/c 2 + 2p • k x + i0 A 

F^ ubl and R\ can be interpreted as the infrared singular and infrared regular parts of the 
auxiliary integrals 

^rd-k^ Kk\ ■ ■ ■ 

(2rr) n [A; 2 - M 2 + zO+]«[A; 2 + 2p-k x + i0+]^ ' 1 ' 

respectively. Since h is a "standard" one- loop integral that is only UV divergent, the addi- 
tional divergences in its IR regular part Ri must cancel exactly with the divergences in its 
IR singular part F 4 ubl > dw , Therefore, 



pSub lt div TDadd 

' , (64) 



e e 
and, using R® dd = —l^ dd ^ it also follows that 

? Sub-i,div T a dd 



rpouo 1 ,aiv j ai 

^ = ^-. (65) 



Analogously 



pSub 2 ,div TDadd jadd 

-± = -il 2 _ = l 2 _. (66) 

e e e 

Having established the relationship between the terms in F 4 and the terms in F 2 and 
F% we now describe the renormalization procedure. Our method consists of treating each 
two-loop integral contributing to F 4 as an independent integral. We then renormalize each 
two-loop integral in the MS scheme, i.e. we 

- determine the divergences in the subintegrals, 

- use the divergences as vertices in one-loop counterterm integrals that are added to F 4 , 
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- perform an additional overall subtraction by absorbing all remaining divergences in 
counterterms, 

- replace jx = ^y/ 2 e 2 ^ - and set jj = m. 

The divergences in the subintegrals are given by F± ubi,div _ The one-loop counterterm 
integrals using these divergences read 

TpSub\,div rpSub2,div 

- ^- F 4 fc2 - ^- F 4 fcl . (67) 

According to Eqs. (1591) . (!60|) . (I65I) . and (|66|) this can be written as 

jadd jadd 

-^h-^h. (68) 



When added to F 4 we obtain 

jadd jadd 

F.-t^-h-tl-h. (69) 

Using the notation of Subsec. IV Al we write F 4 as the sum of the additional divergences and 
a remainder F 4 , 

rpadd.2 T-^add,l 

F 4 = ^V + ^^+F 4 . (70) 
e 2 e 

Note that the divergent terms F^ dd,i are not the divergent expressions F% dd ' 1 of Eq. fl53l) . 
Performing the e expansion for the integrals Ii, 

h = e- 1 lf d + lj° ) + d^ + ---, 
the sum of F 4 and the counterterm integrals is given by 

rpadd,2 rpadd,l j a dd jadd jadd jadd 

b -±— + t± 2^A- - — 4 0) - — 4 0) + F, - lf d I? - lf d I? . (71) 

e 2 e e 2 e e 

We now show that the remaining divergences are analytic in M 2 and can therefore be 
absorbed by counterterms. Recall that the sum of all additional divergences has to vanish, 
since they are not present in the original integral, 



add, 2 j-iadd.l rpadd,2 j-iadd,l rpadd,2 j-iadd,l rpadd,2 rpadd.l 

o = + ^— + ^V- + — + ^V- + — + + — • ( 72 ) 

e 2 e e 2 e e 2 e e 2 



As shown above F 2 and F 3 are the sums of products of one-loop integrals, so Eq. ( 1721) can 
be rewritten as 

j-iadd,2 rpadd,l jadd oadd jadd jyadd jadd jyadd jadd 

o = ^- + ^— + + — 4 0) + A (0) — + + — 4 0) 

e 2 e e 2 e e e 2 e 

jDadd r?add,2 padd,l 

+4 0) ^ + ^V + ^. (73) 
2 e e 2 e 
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Making use of If ,dd = —B% dd the sum of all additional divergences takes the form 

rpadd,2 rpadd,l -padd TDadd 

n = ^ + ^ ^ rV> - *2- rW 



e^ e 

jadd jadd jadd jadd j-jadd,2 j-jadd,l 

2 hA_ _ ^ 7 (o) + + ^ (74) 



All terms in for the two-loop integral as well as the infrared regular terms in one-loop 
integrals are analytic in M 2 . Therefore the first line in Eq. (J71I) is analytic in M 2 . Since the 
sum of all terms vanishes the second line also has to be analytic. This second line, however, 
comprises exactly the remaining divergences in Eq. (ITTj) . which are therefore analytic in M 2 
and can be subtracted. After these divergences have been absorbed in counterterms, the 
resulting expression for the renormalized contribution of F<± reads 

Fl = F 4 -J2 If 1 !? ~ J 2 dd/ i (1) , (75) 



where we have explicitly shown the sums again. Comparing with Eq. ( 1501) we see that our 
result exactly reproduces the expression for the renormalized original integral H\. 



C. e-dependent factors 

For actual calculations it is often convenient to reduce tensorial integrals to scalar integrals 
before performing the dimensional counting analysis as well as the renormalization. The 
reduction of the tensorial integrals can result in e-dependent factors multiplying the scalar 
integrals. These change the form of the result of Eq. (1501) since additional finite terms can 
appear. Let the e-dependent factor be given by 

0( e ) = 0(0 + e 0« + e 2 ^ + ■■■. (76) 

Consider performing the k\ integration first. Suppose that from the result one can extract 
an e-dependent factor (fi(e), and the subsequently performed k 2 integration leads to another 
e-dependent factor, ^(e), with 

0(e) = ^(e)^ 2 (e). (77) 

One can also perform the k 2 integration first, which leads to a different factor <^ 2 (e) , followed 
by the k\ integration resulting in a factor <pi(e) with 

0(e)=^ 2 (e)-^(e). (78) 

The terms v 9 i( e ) = V 9 i°' ) + + e2 fi + ' ' ' an d ^(e) = 02 + e<£>2 + ^fi^ + * ' ' can 
then directly be taken into account when determining the divergent contributions from 
subintegrals. The result for the renormalized integral 0(e) H 2 reads 

hz+ = Ft-^if (vfir + + <pfi?)-<pfir {(pfir + 4 1} 4 0) + ^s o) a (1) ) , 

(79) 

where F% denotes the finite terms in 0(e) F 4 , and (p± and <p 2 ^ are the e-independent 
terms in Jj, <pi and 02, respectively. Our simplified method still holds provided the e- 
dependent factors are taken into account. 
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<b-®-&-(p- 



FIG. 2: Two-loop diagram contributing to the nucleon self-energy. 

As an example consider the diagram of Fig. [2j Ignoring constant factors, one can show 
that in a calculation up to order 0(q 6 ) the nucleon mass only receives contributions from 

-fit ~ + m)Y(p - m)/ J J ^CDE ' (8 ° } 

where the denominators are given in Eq. (1351) . One would also obtain the expression of 
Eq. (j80p if one considered a diagram with fictitious particles as shown in Fig. [3] (a), with 
Feynman rules given by 

j 



k 2 — M 2 + zO" 
i> — m 



P + k k 2 + 2p ■ k + iO H 



f + m 



P + k k 2 + 2p-k + i0+ 

ii 
ii 

T 



7^ 



The subintegral corresponding to performing the k\ integration first is shown in Fig. [3] (b). 
With the Feynman rules above it is proportional to 

(n - 3) (4m 2 7 v - Amp v ) , (81) 

so that we can identify (pi(e) = n — 3 = 1 + 2e. The subsequent k 2 integration corresponds 
to the diagram of Fig. [3] (c), where the diamond-shaped vertex is given by the result of the 
k\ integration. One finds that the term proportional to p u only contributes to higher orders 
and can be ignored. The remaining expression is proportional to 

(n - 3)7^ - mhvgT = -2m(n - 3)(n - 1), (82) 
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(a) (b) 



FIG. 3: Two-loop diagram and diagrams corresponding to k\ and subsequent k 2 integrations. The 
diamond-shaped vertex corresponds to the result of the k\ integration. 

and therefore ^(e) = n — 1 = 3 + 2e. On the other hand, considering the k 2 integration 
first leads to an analogous analysis with the results (p 2 {e) = n — 3 = 1 + 2e = y?i(e) and 
ipi(e) = n- l = 3 + 2e = (p 2 (e). 

In the cases where one cannot identify the individual contributions to 0(e) from the 
integrations of k\ and k 2 , respectively (this happens for example for tensor integrals of the 
type k^k^), one has to perform the dimensional counting analysis before reducing the tensor 
integrals. 

VI. APPLICATION: NUCLEON MASS TO ORDER 0(q 6 ) 

As an application we consider the nucleon mass up to and including order 0(q 6 ). 8 The 



result for the chiral expansion obtained from this calculation has been published in Ref. [19 
Here, we present more details of the calculations. 



A. Lagrangian and power counting 

The effective Lagrangian is given by the sum of a purely mesonic and a one-nucleon part, 
C e ft = C.2 + C± H h C^ N + Cfl + C* N + C^ N + C^ N + Cfl H . (83) 

The purely mesonic Lagrangian at order 0(q 2 ) is given in Ref. Reference 32] contains 
the mesonic Lagrangian at order 0(q 4 ) as well as the lowest-order nucleonic Lagrangian. We 
use the conventions of Ref. [ll| for the Lagrangian at order 0(q 2 ) and of Ref. (34[ for the 
Lagrangians at order 0(q 3 ) and 0(q 4 ). While the complete Lagrangians at order 0(q 5 ) and 
0(q 6 ) have not yet been constructed, up to the order we are considering, vertices from these 
two Lagrangians only appear as contact terms. The light quark masses are proportional 
to the square of the pion mass, and only analytic expressions containing the quark masses 
appear in the effective Lagrangian. Therefore the nucleon mass does not receive any contri- 
butions from the Lagrangian at order 0(q 5 ) in our calculation. The contributions from the 
Lagrangian at order 0(q 6 ) are of the form giM 6 , where g\ denotes a linear combination of 
low-energy coupling constants (LECs) from C n ^. 



Here we consider perfect isospin symmetry. 
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The bare Lagrangians are decomposed into renormalized and counterterm parts. Here 
we only show explicit results obtained from the renormalized Lagrangians, i.e., all appearing 
coupling constants are renormalized coupling constants. The renormalization procedure can 
then be viewed as simply replacing loop integrals by their infrared singular parts. 

We use the following standard power counting [35, Each loop integration in n dimen- 
sions is counted as q n , a pion propagator nucleon propagator as q and vertices 
derived from Li and C^ N as q % and q\ respectively. 



B. Inclusion of contact interaction insertions 

To simplify the calculation we include the self-energy contributions from contact term 



diagrams in the nucleon propagator [111 ]. The advantage of this choice is that all self-energy 
diagrams with contact interaction insertions in the propagator are summed up automatically. 

In terms of the nucleon mass in the chiral limit m the full nucleon propagator can be 
written as 

S N {f) = * ... n+ , (84) 

where — i£(i), m) is the sum of all one-particle irreducible self-energy diagrams. The physical 
nucleon mass is determined by the solution to the equation 

s N\ =mN = \t-™-ntM\\ f=mN = v. (85) 

The self-energy receives contributions from contact terms as well as from loop diagrams, 

£Qz>, m) = S c + Zioop^p, m). (86) 

Due to the form of the BChPT Lagrangian used here, E c for the nucleon is independent of 
ti>. Inserting Eq. flHBI) into Eq. fl85|) one finds 

\f - m - E c - T, loop (f, m)] \ f=mN = 0. (87) 

In a loop expansion Eq. (I8"7j) has the perturbative solution 

m N = m + J: c + 0(h). (88) 

In the above the propagator which is used in the calculation of the self-energy diagrams has 
been chosen to be 

1 i> — m + iu + 
However, one can also choose this propagator to be 

^(p) = ~ V - - n+ ■ (90) 

This corresponds to including in the free Lagrangian those terms bilinear in \l/ which 
generate the contact term diagrams in the self-energy contribution. The advantage of this 
choice is that all self-energy diagrams with contact interaction insertions in the propagator 
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are summed up automatically. With this choice of the propagator the self-energy is now 
given by the sum of loop diagrams only, i.e. 

m) -> Eioopif, fh) , (91) 

where 

fh = m + S c . (92) 

As an additional benefit, when working to two-loop accuracy, one can set i> = fh in the 
expression of two-loop diagrams, since corrections are at least of order 0(h 3 ). To obtain 
m N one has to solve the equation 



=mjv 



S N l (m N ) = \jf>-fh- Sjoop (^,m)] ^_ =0. (93) 
Inserting the loop expansion for E ioop (^, fh) 



Z loop (± fh) = Jffifi-fo fh) + h 2 nfU± fh) + ■ ■ ■ , (94) 



using the ansatz 

mjy = fh + hAmi + h 2 Arri2 + • • • (95) 
and expanding around fh we obtain up to the two-loop level 



= fh + hAm 1 + h 2 Am 2 — m — fiE^(m + Mmi,m) — /i 2 s|^ (m, m) 



+ £ 2 



Am 2 - Ami S^(m, m) - S^ p (m, fh) 



(96) 



where s|^'(|l, m) denotes the derivative of S^ {jf>, fh) with respect to f. The solutions for 
Ami and Am 2 are given by 

Am i = s£Jp(m,m), (97) 
Am 2 = E£Jp(m, m)S[^(m, m) + s[ 2 ] p (m, m) . (98) 

To obtain the nucleon mass up to chiral order 0(q 6 ) one needs to determine S c , Ami 
and Am 2 up to that order. In the following we will not directly evaluate Am 2 as given 
in Eq. fl98|) since we are only infested in the combination hAmi + /i 2 Am 2 with h — 1. 
Instead, as indicated in the first line of Eq. fl96l) . we will use the result for Am x to determine 
S|^ p (m + hAm 1} fh) directly (see Eq. H103[) below). 

In principle, the nucleon propagator is a 2 x 2-matrix in isospin space. For arbitrary 
values of the up and down quark masses the propagator is a diagonal matrix; since however 
in this work the isospin-symmetric case m u = rrid is considered, the masses of proton and 
neutron are identical and the propagator is proportional to the unit matrix. 



C. Contact terms 

The contributions to the nucleon mass from contact interactions are given by 

5m c = -4ciM 2 - (16e 38 + 2e U5 + 2en 6 )M 4 + g x M % 

= -4ciM 2 - eiM 4 + 3iM 6 , (99) 
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<D b <b — + — CD (£> (6 G> 

(a) (b) (c) 




(d) (e) (f) 

FIG. 4: One-loop diagrams contributing to the nucleon self-energy up to order 0(g 6 ). 

where M 2 is the lowest-order expression for the square of the pion mass. We use the notation 

e x = 16e 38 + 2eii 5 + 2e 116 (100) 
and gi denotes a linear combination of LECs from the Lagrangian at order 0(q 6 ). 

D. One-loop diagrams 

The one-loop diagrams contributing to the nucleon mass up to order 0(q 6 ) are shown in 
Fig. (j3J). Diagrams (a) and (d) are of order 0(q 3 ) and 0(q 4 ), respectively, and have been 
determined in Ref. Diagrams (b) and (c) are of order 0(q 5 ), while the power counting 
gives D = 6 for diagrams (e) and (f). 

Using dimensional regularization the unrenormalized results for the one-loop diagrams 
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up to order 0(q 6 ) read 



Si(fe) 



3gi 
4F 2 



+ m)#iv + (p + m)M*H„ N (p 2 ) + (p 2 - m 2 )fH^(p 



r(p) 



3g 



.4 



(2d 



_ 16 -d 18 )M 2 
+(p 2 - fh 2 )fH^{p 2 ) 



m)M 2 H wN (p 2 ) 



2F 4 



+ m)M 2 / 3 M 2 - (p 2 - m 2 )H^ N (p 



TiV 



H 1T -{p -m )H nN (p 



+ m)H. 



(p) 

txN 



3 

F 2 " 



(2 Cl -c 3 )M 2 H n -c 2 ^H^ 



{ [2(e 14 + e 19 ) - e 36 - 4e 38 ] M 4 ^ + 2 [e 15 + e 20 + e 35 ] ^M 2 H^ 



+6e 16 4^r oo) i 



S K/) = 6 ^4 \ 2 Cl [l 3 M 2 H^ - hH w ] ~-— 2 [{h ~ h)H n + / 3 M 2 #_] 

-c 3 [(/ 3 - / 4 )^ - Z 3 M 2 #_] } . (101) 

The integrals H n , H nN (p 2 ), . . . are given in App. [Bl Various combinations of fourth-order 
baryonic LECs appear through the vertex in diagram (e). To simplify the notation we use 
ei as defined in Eq. (11001) and 



e 2 = 2eu + 2ei 9 - e 36 - 4e 38 
£3 = ei 5 + e 20 + e 35 



(102) 



for these combinations in the following. 

To determine the contribution of these diagrams to the nucleon mass we evaluate the 
expressions of Eq. (11011) at p 2 = m 2 N between on-shell spinors. To the order we are working, 
we can use = fh + hArrii and thus obtain £[*^(m + hAmi,rh) in Eq. (1961) . Further, 
we renormalize the one-loop integrals by replacing them with the corresponding infrared 
singular parts. The infrared renormalized expressions for the mass contributions, denoted 
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by a superscript r, up to order M 6 are given by 
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The scale dependence of the renormalized low-energy constants is governed by 

k,o = k(n) + Ji\ d ifl = di(fi) + -^A, e i:0 = ei(fi) + —^2 A, 
where the subscript denotes bare quantities and 

A = £{^2> l(4 " + r ' (1 > + 4 
The coefficients 7$ are given by 

1 

7s = 74 = 2, 



while the & can be taken from Ref. 33 



Note that the numbering of terms in the Lagrangian used here differs from Ref. [3 
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The LECs <32,e3,ei6 appear through the vertex shown in Fig. @] (e), which also gives the 
contact term contribution at order 0{q i ) to ttN scattering as analyzed in Ref. j37|. We can 
therefore relate the combinations of LECs used here to the ones in Ref. 37j (here denoted 
by a superscript BL), resulting in 



ei 
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"1 5 



Using the expressions for the renormalized couplings given in App. E of Ref. [37| we find 
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1 + 3g^ + yg^ + 8cim + c 2 m - 4c 3 m 



52 

10 + 12g^ + — g\ + 8c 2 m 



(107) 



E. Two-loop diagrams 

The two-loop diagrams relevant for a calculation of the nucleon self-energy up to order 
0(q 6 ) are shown in Fig. [5j According to the power counting there are further diagrams at 
the given order. An example would be diagram [5] (c) with one first-order vertex replaced by 
a second-order one. As a result of our calculation we find that these diagrams give vanishing 
contributions to the nucleon mass up to the order we are considering. 

We again employ dimensional regularization. The unrenormalized expressions for the 
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FIG. 5: Two-loop diagrams contributing to the nucleon self-energy up to order 0{q % ). 
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mass contributions of the diagrams of Fig. [5] up to order 0(q 6 ) can be reduced to 
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The integrals H^^ia, b\n) and H 2 (a, b, c, d, e\n) are defined in App. O Here we have ex- 
pressed tensor integrals in terms of scalar integrals in higher dimensions where convenient 
(see App. IO and also Ref. (38|). 

After performing the infrared renormalization as described in Sec. |V]the contributions to 
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the nucleoli mass up to order 0(q e ) read 
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F. Results and discussion 



Combining the contributions from the contact interactions with the one- and two-loop 
results we obtain for the nucleon mass up to order 0(q 6 ) 

m N = m + k x M 2 + k 2 M 3 + k 3 M 4 In — + k 4 M A + k 5 M 5 In — + A; 6 M 5 

A* A* 

+ k 7 M 6 In 2 — + k 8 M 6 In — + A; 9 M 6 . (110) 
// A* 



29 



The coefficients ki are given by 
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In general, the expressions of the coefficients in the chiral expansion of a physical quantity 
differ in various renormalization schemes, since analytic contributions can be absorbed by 
redefining LECs. However, this is not possible for the leading nonanalytic terms, which 
therefore have to agree in all renormalization schemes. Comparing our result with the 
HBChPT calculation of 2(j, we see that the expressions for the coefficients fc 2 , fc 3 , and 
k 5 agree as expected. At order 0(q 6 ) also the coefficient k 7 has to be the same in all 
renormalization schemes. Note that, while k§M^ and k$M 6 In — are nonanalytic in the quark 
masses, the algebraic form of the coefficients he and k$ are renormalization scheme dependent. 
This is due to the different treatment of one-loop diagrams in different renormalization 
schemes. The counterterms for one-loop subdiagrams depend on the renormalization scheme 
and produce nonanalytic terms proportional to M 5 and M 6 In — when used as vertices in 
counterterm diagrams. We find that our result for k§ coincides with the HBChPT calculation 



of Ref. [20] except for a term proportional to <i 2 g, which, however, does not have a finite 
contribution for manifestly Lorentz-invariant renormalization schemes [39]. Therefore, at 
order 0(q 5 ) the chiral expansion of the IR renormalized result reproduces the HBChPT 
result. 

The result for the nucleon mass should be scale-independent at each order in the chiral 
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expansion, and showing this scale independence serves as a check of our results. The terms 
up to and including order 0(M A ) have been discussed previously (37| . Using the expressions 
for the scale dependence of the renormalized couplings of Eqs. (11041) . (I105j) and f II 6 j) we 
see that the contribution at order C(M 5 ) is independent of \i as required. At order C(M 6 ), 
the scale dependence of the LEC g\ is not known which prevents a complete analysis at 
this order. However, gi does not contribute to terms proportional to In M In /i since it 
only appears in the analytic expression at 0(M 6 ). We can therefore analyze the terms 
proportional to In M In fi which must vanish if our result is to be scale- independent. This is 
the case, which can be shown using the expressions of Eqs. (I104p - (I107I) . 

The numerical contributions from higher-order terms cannot be calculated so far, since 
most expressions in Eq. (11 111) contain LECs which are not reliably known in IR renormal- 
ization. In order to get an estimate of these contributions we consider several terms for 
which the LECs have previously been determined. The coefficient k 5 is free of higher-order 
LECs and is given in terms of the axial-vector coupling constant g A and the pion decay 
constant F. While the values for both and F should be taken in the chiral limit, we 
evaluate k 5 using the physical values g A = 1.2695(29) [40] and F w = 92.42(26) MeV. Setting 
H = m N , m N = (m p + m n )/2 = 938.92 MeV, and M = M n + = 139.57 MeV we obtain 
k$M 5 \rt(M /m^) = —4.8 MeV. This amounts to approximately 31% of the leading nonan- 
alytic contribution at one-loop order, /^M 3 . The mesonic LECs appearing in k§ can be 
found in Ref. [2( and are given by I^itln) = 1.4 x 10~ 3 and I^itln) = 3.7 x 10" 3 at the scale 
/i = m N . The parameter <i 18 can be related to the Goldberger-Treiman discrepancy [37| and 
is given by di% = —0.80 GeV" 2 . The LEC d\Q, however, is not as reliably determined. In or- 
der to estimate the magnitude of the contribution stemming from fc 6 we use the central value 



from the reaction irN — > irirN, di 6 (m N ) = — 1.93GeV |4l|, |42j. It should be noted that 



the calculation of Ref. [41[ was performed in HBChPT, and employing the obtained value 
for in an infrared renormalized expression therefore only gives an estimate of the size of 
the corresponding term. The resulting contribution is kgM 5 = 3.7 MeV and cancels large 
parts of the nonanalytic term k$M 5 ln(M/mjv). In Ref. [43j the parameter d\Q has been de- 
termined by a fit to lattice data. At the scale [i = uin it is given by d^m^) = 4.11 GeV" 2 , 
which does not agree with the result from the reaction tcN — > tttcN. With this value oidr\ 6 
we find fc 6 M 5 = —7.6 MeV. The LECs appearing in k 7 have been determined in Ref. 371 ] . 
and we obtain k 7 M 6 In 2 (M/m N ) = 0.3 MeV. 

The terms k s and kg contain LECs from the fourth order Lagrangian C^ N which have 
not been determined. We try to get a very rough estimate of the size of these contributions 
by assuming that all these LECs as well as g\ are of natural size, that means ~ 1 GeV~ 3 
and gi ~ 1 GeV~ 5 . We choose the d±Q value from nN — > nnN and use the above values for 
the other LECs. Setting all appearing = GeV -3 gives a contribution k$M 6 ln(M J itln) w 
10~ 2 MeV. The choice e { = 5 GeV" 3 results in k 8 M 6 \n(M/m N ) w 0.9 MeV, while e { = 
— 5 GeV~ 3 gives k$M 6 ln(M/mjv) ~ —0.9 MeV. A similar analysis for the term kgM 6 gives 
k 9 M 6 w -2.8 MeV for all ej = GeV" 3 and g x = 0GeV~ 5 , while setting e { = 5 GeV" 3 , 
gi = 5 GeV" 5 and e { = -5 GeV" 3 , g x = -5 GeV" 5 results in k g M 6 ss -2.5 MeV and 
fc 9 M 6 w —3.2 MeV, respectively. One should note, however, that the numbers obtained 
here are only very rough estimates. Choosing eu = ei 5 = 5 GeV -3 and ei 6 = eig = e 20 = 
635 = = = 1 GeV -3 , g\ = 1 GeV -5 leads to large cancelations between the terms 
fcgM 6 ln(M/m.jv) and kgM 6 , resulting in a complete contribution at order 0(M 6 ) of about 
0.3 MeV. As a check we also use the value of die a s obtained in Ref. (43[, which results in 
contributions from kg that are about a factor 10 larger, while the dependence of kg on d\% 
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FIG. 6: Pion mass dependence of the term k$M 5 ln(M/mjv) 
comparison also the term k 2 M 3 (dashed line) is shown. 



(solid line) for M < 400 MeV. For 



is much less pronounced. Clearly a more reliable determination of the higher-order LECs is 
desirable. 

Chiral expansions like Eq. (11 101) play an important role in the extrapolation of lattice 
QCD results to physical quark masses, and the nucleon mass is an example that has been 



studied in detail (see, e.g., Refs. (H, a, EBI, |47j ) . In Ref. |44| su ch an extrapolation was 
performed for the nucleon mass up to order 0(q 4 ), while Ref. "0 includes an analysis of 



the fifth-order terms. It was shown, as had also been argued in Ref. |42(], that the terms 
at order 0(q 5 ) play an important role in the chiral extrapolation. As an illustration we 
consider the leading nonanalytic term at this order, k^M 5 ln(M/mjv)- Its dependence on 
the pion mass is shown in Fig. [6] for pion masses below 400 MeV, which is considered a 
region where chiral extrapolations are valid (see, e.g., Refs. (48|, H^]). We see that already at 
M ~ 360 MeV the term k$M b hi(M / tun) becomes as large as the leading nonanalytic term 
at one-loop order, k 2 M 3 , indicating the importance of the fifth-order terms at unphysical 
pion masses. Since the contribution at order 0(M 6 ) depends on a number of unknown LECs, 
we do not attempt to perform a chiral extrapolation up to this order here, but restrict the 
discussion on the pion mass dependence of the term kjM % \v? (M / m^) ■ Figure [7] shows 
this dependence for pion masses below 400 MeV. No errors are given for the LECs Ci, c 2 , 
and C3 in Ref. 
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For an estimate we have assumed the relative errors of these LECs 
and of g A to be 20%, and the corresponding error for k 7 M G ln 2 (M/mjv) is shown in Fig. [71 
For comparison we also show the nonanalytic term at fourth order, /C3M 4 \n{M/m^). As 
expected, and in contrast to the fifth-order term, the two-loop term kjM & In 2 (M/mjy) is 
smaller than the one-loop contribution fc 3 M 4 ln(M/m^) in the considered pion mass region. 
Note that the relative difference in the pion mass dependence between fc 5 M 5 \n(M/m^) 
and k 2 M 3 , as well as kjM 6 In 2 (M/m^v) and k 3 M 4 hi{M/m^) is proportional to a factor 
M 2 \n(M/m^f), and that for the physical pion mass the differences in the two cases are 
comparable on an absolute scale. We also show the pion mass dependence of the terms 
k 7 M 6 In 2 (M/m N ) and k 3 M A ln(M/m N ) up to M « 700 MeV, which, however, is beyond the 
domain that is considered suitable for the application of Eq. (II 101) . Again the sixth-order 
term remains much smaller than the fourth-order one, also at higher pion masses. However, 
the above considerations are not reliable predictions for the behavior of the complete two- 
loop contributions at unphysical quark masses. This is because here only one of the terms 
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FIG. 7: Pion mass dependence of the term k-jM® In 2 (M/rajy) (solid line) for M < 400 MeV. The 
shaded band corresponds to relative errors of 20% in the LECs. For comparison also the term 
k^M 4 ln(M/mj\r) (dashed line) is shown. 
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FIG. 8: Pion mass dependence of the term k-jM® lii 2 (M/m,]\f) (solid line) for M < 700 MeV. The 
shaded band corresponds to relative errors of 20% in the LECs. For comparison also the term 
A^M 4 m(M/mjv) (dashed line) is shown. 

at order 0(q 6 ) is considered, and the contribution of the analytic term proportional to M 6 
can be considerably larger than k 7 M 6 \v?(M /m^) depending on the values of the unknown 
LECs. 



G. Nucleon a term 



The Feynman-Hellmann theorem [5TJ, [51 



relates the nucleon mass to the value of the 
nucleon scalar form factor at zero momentum transfer, the so-called a term (see, e.g., [52|, 



a(q 2 = 0) = M 



,drriN 
dM 2 



(112) 
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Applying the Feynman-Hellmann theorem to Eq. (11 101) . the chiral expansion of cr(0) is given 

by 

<r(0) = hM 2 + | k 2 M 3 + 2k 3 M A + 2kA M 4 + - k 5 M 5 In — 

+ - (k 5 + 5k 6 )M 5 + 3k 7 M 6 In 2 — + (k 7 + 3k 8 )M 6 In — + ( — + 3k 9 ] M 6 . 

(113) 



The first four terms have already been determined in Ref. [llf . To estimate the contributions 
of the terms of order 0(M 5 ) we use the same values for the LECs as above, in particular 
the value of die as extracted from the reaction ttN — > ttttN. The combined contributions at 
order 0(M 5 ) are 

- k 5 M 5 In — + - (fc 5 + 5k 6 )M 5 « -0.2 MeV. (114) 
2 1 1 2 

Due to the dependence of the order 0(M 6 ) nucleon mass contribution on the specific values 
of the unknown LECs e^, we do not attempt to evaluate the terms at order 0(M 6 ) in 
Eq. (PS). 



VII. SUMMARY 



We have shown details of how to consistently renormalize two-loop diagrams in manifestly 
Lorentz-invariant BChPT within the framework of infrared regularization. The renormal- 
ization procedure preserves all relevant symmetries such that renormalized expressions fulfill 
the relevant Ward identities. Renormalized diagrams also obey the standard power count- 
ing. We have presented a simplified method of renormalizing diagrams with one small scale, 
which relies on dimensional analysis. Integrals of this kind appear, e.g., in the calculation 
of the nucleon mass or the axial- vector coupling constant. In this method integrals derived 
from the original expressions are renormalized using the MS scheme, which simplifies the 
calculations considerably. As an application we have calculated the nucleon mass up to and 
including order 0(q 6 ). For physical values of the pion mass, the numerical estimate of the 
two-loop contributions is reasonably small. For example, the estimate of the 0(q 5 ) term of 
the nucleon a term is —0.2 MeV. However, when considering the nucleon mass as a function 
of the pion mass, we have seen that already at a pion mass of 360 MeV the nonanalytic 
contribution at order 0(q 5 ) may become as large as the nonanalytic 0(q 3 ) contribution. 
From this one cannot conclude that the chiral expansion breaks down at this value of the 
pion mass since the analytic terms at 0(q 5 ) might cancel parts of the nonanalytic term. One 
should, however, take special care when performing chiral extrapolations beyond this value. 
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APPENDIX A: DIMENSIONAL COUNTING ANALYSIS 



Analytic expressions for two-loop integrals, especially when two mass scales such as the 
pion mass M and the nucleon mass in the chiral limit m appear in the same integral, can be 
extremely difficult to obtain. Since we are interested in the chiral expansion of the considered 
integrals in the present work, we do not have to find a closed-form solution to the appearing 
integrals, but can use a method called dimensional counting analysis [30[ for the evaluation 
of integrals. A closely related way of calculating loop integrals is the so-called "strategy of 



regions" [54(. Here we present an illustration of dimensional counting for one- and two- loop 
integrals. 



1. One-loop integrals 

The advantage of the dimensional counting analysis for one-loop integrals lies in its ap- 
plicability to dimensionally regulated integrals containing several different masses. Consider 
integrals with two different mass scales, M and m, where M < m, and a possible external 
momentum p with p 2 « m 2 . Dimensional counting provides a method to reproduce the 
expansion of the integral for small values of M at fixed p 2 — m 2 . To that end one rescales 
the loop momentum k i— > M ai k, where a, is a non- negative real number. After extracting an 
overall factor of M one expands the integrand in positive powers of M and interchanges sum- 
mation and integration. The sum of all possible rescalings with subsequent expansions with 
nontrivial coefficients then reproduces the expansion of the result of the original integral. 

To be specific, consider the integral 

H ^ip ) = J ( fc 2 _ m 2 + iO+)[(k + p) 2 -m 2 + i0+] ' ^ 

It can be evaluated analytically and the result is given in App.[B] After rescaling one obtains 

j 0\ if M na *d n ~k IK , 

H nN (p 2 ) i — ► - — — / -= = = . A2 

(2tt)« J [k^M 2 ^ — M 2 + ^0+] [k 2 M 2a i + 2p ■ kM^ + p 2 - m 2 + i0+] 

No overall factor of M can be extracted from the second propagator, which is therefore 
expanded in positive powers of M. As a result the integration variable k only appears in 
positive powers in the expanded expression of this propagator. If < ai < 1 one can extract 
the factor M~ 2oLi from the first propagator, which takes the form 

(A3) 



k 2 - M 2 ~ 2 ^ + i0+ 



Expanding in positive powers of M and interchanging summation and integration one obtains 
integrals of the type 

d n k^ — . (A4) 

(k 2 + io+y 

Combined with the expansion of the second propagator the resulting coefficients in the 
expansion in M are integrals of the type 

k m 

d n k^— , (A5) 

(k 2 + io+y 
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which vanish in dimensional regularization. For the case 1 < a, the first propagator in 
Eq. (1A2j) can be rewritten as 

' ' (A6) 



M 2 (k 2 M 2a ^ 2 - 1 + i0+) 



Expanding in M and combining with the expansion of the second propagator one obtains 
integrals of the type 

r d n k k j , (A7) 



which, again, vanish in dimensional regularization. The only contributions to H n N(p 2 ) can 
therefore stem from a,- = and a,- = 1. For a,- = one obtains 



**nKP ) (2tt)-^ 1 > J [k 2 + tO+Y+^k + p) 2 - m 2 + 10+] 1 { ' 

while the expression for 014 = 1 reads 

F d) ,*s _ _±_ fv lY M n ~ 2+i r <rHk?M + 2 P .~ k y 

Kn[P } ~ (2k)- ^ l) (P 2 - m 2 )^ J p - 1 + *0+] ■ { } 

The expansion of H 1t n(p 2 ) is then given by 

hMp 2 ) = h^(p 2 ) + h^( p 2 ), (A10) 

which correctly reproduces the result of App. [Bl 



2. Two- loop integrals 

While one of the advantages of the dimensional counting method lies in its applicability to 
integrals containing several mass scales, a difficulty arises for the calculation of the nucleon 
mass. Since integrals have to be evaluated on-mass-shell, the two small scales M and p 2 — m 2 
are not independent of each other and are comparable in size. Therefore an expansion in 
2^f m 2 does not converge. By the choice of the nucleon propagator mass to include all 
contact interaction contributions, the terms p 2 — m 2 in the propagator can be neglected 
in two-loop integrals since they are of higher order in the loop expansion. The two-loop 
integrals contributing to the nucleon mass are therefore reduced to integrals with only one 
small mass scale, for which an expansion in M can be obtained. 

For the extension of the dimensional counting method to two-loop integrals 



H 2 (a, b, c, d, e\n) 



d n k x d n k 2 1 



' ' ' ' 1 ' J J (2tt) 2 ™ [k 2 - M 2 + i0+] a [k 2 - M 2 + i0+] b 

1 

pfcf + 2p ■ ki + iO+] c m + 2p-k 2 + iO+YKh + k 2 ) 2 + 2p ■ k x + 2p ■ k 2 + iO+Y ' 

(An: 
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one has to consider all possible combinations of rescaling the integration variables k\ \— > 
M ai ki, k 2 i— >• M^k 2 . The expansion of the two- loop integral is then given by 

H 2 (a, b, c, d, e\n) = M^ aiA) h {aiA \a, b, c, d, e\n) 

= J2 H {a * A \a, b, c, d, e\n), (A12) 

where ip(a>i, fy) is the overall power of M extracted for each rescaling, the functions 
h^ A \a,b, c,d, e\n) are the expressions for the integrated expansions, and we have defined 
H^ aul3i \a, b, c, d, e\n) = M f, ' ai ' ft ^' a,,ft '(a, b, c, d, e\n) to simplify the notation. Following the 
discussion of the one-loop sector one sees that the only combinations that give non- 

vanishing contributions are (0, 0), (1, 0), (0, 1) and (1, 1), so that a two-loop integral is given 
by 



H 2 (a, 6, c, d, e\n) = H {0 ' 0) (a, b, c, d, e\n) + H^'°\a, b, c, d, e\n) + H {0 ' 1] (a, 6, c, d, e 



n 



+ # (1>1) (a,6,c,d,e|n). (A13) 

To shorten the notation and to avoid confusion with other superscripts used for further 
expansions in this paper the contributions corresponding to H^°'°\ H^'°\ H^ ' 1 ', and H^ 1 ' 1 ^ 
are also denoted by F\, F 2 , F 3 and F 4 , respectively 

From a technical point of view it is convenient to consider the rescaling k± h- > (M/m) ai ki, 
since then the integration variables k have dimension of momenta. This 
also facilitates the evaluation of certain loop integrals appearing in the calculation of the 
nucleon mass. 

As an example consider the integral H 2 (l, 1, 1, 1, l\n). For (0,0) the resulting integrals 
read 



x 



d n k x d n k 2 1 



(2tt) 2n [k\ + iQ+} l+i [kl + i0 + } l+ i[kf + 2p ■ k x + 
1 



[k\ + 2p-k 2 + z0+] [(h + k 2 ) 2 + 2p-k 1 + 2p-k 2 + i0+] ' 

(A14) 

While still a two-loop integral that does not directly factorize into the product of one-loop 
integrals, the vanishing of the mass scale M simplifies the evaluation of the integral. The 
rescaling of only ki leads to 



3+2i+j+l irrl 3-n-j-l 



d n k x d n k 2 1 



(k!y& 1 + ^p-h + 2k 1 -k 2 y 



\k\ + i0+] 1+i [2p ■ ki + tO+Y+im + 2p-k 2 + i0+] 2 + l 

(A15) 

while the expression for H^Hl, 1,1, 1,11 n) can be obtained by substituting ki i— > k 2 and 
k 2 I— > ki in Eq. (IA15|) . One sees that the integrals of Eq. (IA15I) can be reduced to the 
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product of tensorial one-loop integrals, which is a considerable simplification compared to 
the original integral. The last contribution stems from = 1, /?» = 1 and reads 



M 

<-./.' 



d n k x d n k 2 



1 



X — 



[A; 2 - m 2 + i0 + ][2p ■ ki + i0+] 1+i [2p ■ k 2 + i0 + ] 1+ i[2p ■k 1 + 2p-k 2 + i0+] 1+l ' 

(A16) 

where the integration can be reduced to the evaluation of a set of basis integrals (see 
App. ID]) . The sum of all four contributions reproduces the M expansion of the integral 
F 2 (l,l,l,l,l|n), 

H 2 (l, 1, 1, 1, l|n) = H^°\l, 1, 1, 1, l|n) + ff^l, 1, 1, 1, l|n) + ^(l, 1, 1, 1, l|n) 

+ if (1 ' 1) (l,l,l,l,l|n). (A17) 



APPENDIX B: INTEGRALS AT THE ONE-LOOP LEVEL 



Using dimensional regularization [55( the one-loop integrals are defined as 

H w = 



pu pa I up va , n p& n vp 
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grg +grsr 
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pv rr(00) 
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k^kPk* 
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(27r) n 
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[k 2 
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(27r) n 
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-M 2 + z0+] 2 ' 
1 




(2n) n 
d n k 


k 2 - 


- m? + z0+ ' 
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(2?r) n 
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[k 2 


- M 2 + zO + ][(fc + p) 2 - 
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(2?r) n 


[k 2 


- M 2 + i0+}[(k + p) 2 - 
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- m 2 + ?0+] ' 


(27r) n 
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The tensorial loop integrals can be reduced to scalar ones 56[ and we obtain 

M 2 



#(oo) 



n 



#(0000) 



M 4 
n(n + 2) 



2p 2 
1 



[H v -H N - (p 2 -m 2 + M 2 )H« N (p< 



[H nn - H nN - (p 2 - m 2 + M 2 )H^ N (p< 



Defining 



A 



n— 4 



1 1 



f 6vr 2 n - 4 2 



[ln(4vr) - lE + 1] 



where 7^ = — r'(l) is Euler's constant, and 

p 2 - m 2 - M 2 



n 

the scalar loop integrals are given by [l5| 

M 2 , M 

rr ln— ' 



2mM 



H n = 2M Z A + 



H W7r = 2A + 



16tt 2 



l + 21n^' 



Fjv = 2m 2 A + 



9 , m , m 
~ x -In—. 



57T' /i 



H vN {p 2 ) = 2A + 



16tt 2 



p 2 — m 2 + M 2 , M 2mM ln ' 
-1 + ^ In — + — —F(Q) 

p A jj, p A 



where 



yo 2 - 1 In (-n - v^ 2 - 1) , < -1, 

= { VI -P 2 arccos(-fi), -1 < O < 1, 

V^ 2 - 1 in (n + vn 2 - 1) - w^ 2 - 1, 1 < fi. 

The integral _/v can be obtained from H nN (p 2 ) by differentiating with respect to M 2 



APPENDIX C: INTEGRALS AT THE TWO-LOOP LEVEL 



We define the one-loop integral 



H vN (a,b\n) 



d n k 



(2ir) n [k 2 -M 2 + i0+] a [k 2 + 2p ■ k + i0~ 



39 



Note that we have not included a factor i in the definition, since the one-loop integrals 
HnN{a>-,b\n) always appear in the product H^Nicti, bi\ni)H n jsr(a 2 , b 2 \n 2 ). 
The two-loop integrals #2(0, b, c, d, e\n) are defined as 



H 2 (a, b, c, d, e\n) 



d n k 1 d n k 2 1 



(27r) 2 ™ A a B b C c D d E e 



where 



A = k 2 -M 2 + iO + , 
B = kl-M 2 + iO + , 
C = k 2 + 2p-k 1 +i0 + , 
D = k\ + 2p- k 2 + i0 + , 

E = k 2 l +2p-k 1 + 2k 1 -k 2 + 2p-k 2 + k 2 2 + i0 + . 

The product of two one-loop integrals with the same space-time dimension n can then be 
written as 

H^ N (ai, b 1 \n)H^ N (a 2 , b 2 \n) = H 2 (a u a 2 , 61, b 2 , 0|n). 

We do not attempt to evaluate the integrals H 2 (a,b,c,d,e\n) here. Instead the expres- 
sions for the integrals relevant to the nucleon mass including the corresponding counterterm 
contributions are given in App. ID1 

Tensorial integrals have been reduced to scalar ones in the same dimension using methods 
similar to the one-loop integrals, or the following relations to scalar integrals in higher 
dimensions have been used: 



H^i^a, b, c, d, e\n) = J J 
= -16vrV 



d n k x d n k 2 kf 



(2n) 2n A a B b C c D d E e 
bcH 2 (a, b + 1, c + 1, d, e\n + 2) + cdH 2 (a, b, c + 1, d+ 1, e|n + 2) 

+ ce H 2 (a, b, c + 1, d, e + l\n + 2) + be H 2 (a, b + 1, c, d, e + l\n + 2) 
d n hd n k 2 k^k v x 



H^'fa, b,c,d,e\n 
(An) 2 



(2n) 2n A a B b C c D d E e 
[bH 2 (a,b + l,c,d,e\n + 2) + dH 2 (a,b, c, d + 1, e\n + 2) 



2 

+e H 2 (a, b, c, d, e + l\n + 2)] + <D(p), 

d n hd n k 2 kfk»k$ 



H% uX '(a, b,c,d,e\n 



(2vr) 2 ™ A a B b C c D d E e 



{ '" ' 1 [g^p x + g^p v + <fV] [6 (6 + 1) cH 2 (a, b + 2,c + l,d, e\n + 4) 



2 

+ 6 (6 + 1) e # 2 (a, b + 2, c, d, e + l|n + 4)] + C(p 3 ), 

ff d n hd n k 2 gfc| 
^ (a,b,c,d,e\n) = JJ ^ ^^E* 

= -^r- 9^ e H 2 (a, 6, c, d, e + l|n + 2) + 0(p), 
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H 2 ' a ^{a, b, c, d, e\n) — J J 



d n k x d n k 2 Kk%k p 2 



( 4?f ) -,aPu 



2 



(2tt) 2 ™ A a B b C c D d E e 
g ali p^ [cH 2 (a, b,c+l,d,e\n + 2) + eH 2 (a,b } c, d, e + l\n + 2)] 



+ [<7 a V + 0"V + 9^P a ] [ae(e + 1) H 2 (a + 1, b, c,d,e + 2\n + 4) 

+ adeH 2 (a + 1, b, c, d + 1, e + l|n + 4) 
+ cde H 2 (a,b,c + l,d+ 1, e + l|n + 4) 
+de(e + l)H 2 (a, b, c, d + 1, e + 2|n + 4)] + 0(p 3 ), 

d n k 1 d n k 2 k^k{k^k v 2 



H2^ ,fiu (a,b,c,d,e\n) = J J 



( 2vr )2n A a B b C c D d E e 



= [9 aP 9" u + sTg^ + 9 av 9^\ e (e + 1) H 2 (a, b,c,d,e + 2\n + 4) 

+ g^g^H^a, 6, c, d, e\n + 2) + 0(p). 

Here, 0(p) stands for terms proportional to p p , where p denotes the Lorentz index corre- 
sponding to the integral under consideration. In our calculation of the nucleon mass these 
terms appear in combination with expressions like (p — m) , y p (p+m), resulting in higher-order 
contributions that are not considered. 



APPENDIX D: EVALUATION OF TWO-LOOP INTEGRALS 

As seen in Sec. |V]the calculation of the two-loop integrals relevant to the nucleon mass 
reduces to the evaluation of the F 4 part of the respective integrals. The F 4 parts are sums 
of tensor integrals, which can be reduced to scalar integrals 5|| of the form 



H2~ ,1 \a, b, c, d, e\n)— 1 ^ 1 



(2tt) 2 " [k\ - m 2 + iO+] a m - m 2 + i0 + ] b [2p -kt+iO 



+ c 



X [2p ■ k 2 + i0+] d [2p ■k 1 + 2p-k 2 + i0+} e ' ^ 

where the superscript (1,1) indicates that these integrals have been obtained after rescaling 
both integration variables (see App. [A} and a, b, c, d, e are integers. Depending on the values 
of the exponents c, d, and e, one can evaluate H^' 1 \a, b, c, d, e\n) with the help of several 
basic integrals. 



1. e = 

If the exponent e vanishes, the integral can be written as the product of one-loop integrals, 
H£ A) (a, b,c,d,0\n) = H^\a lC \n)H^ ] {b,d\n), (D2) 
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where 



H[ 1] (a,b\n) 



d n k 



i 



(2tt)™ [k 2 - m 2 + i0+} a [2p ■ k + iO+] b 

l-2o-26 



r[|]r[a + |-f] 



2 6( 47r )n/ 2 r[a]r[^±i] 



[m 



2^n/2-a-b/2 ^2\-6/2 



(D3) 



2. c = d = 0,e/0 

If c = d = 0, the expression for the integral H 2 (a, b, 0, 0, e\n) reads 

H£ ,l \a,b,Q,0,e\n) 

d n hd n k 2 1 

(27r) 2n \k\ -m 2 + i0+} a [k 2 - m 2 + i0+] fe [2p ■ ^ + 2p ■ fc 2 + iO + Y 



;2-2a-2b-2e 



2(4vr) 



( TO 2)n-a-6-e/2 (p 2)-e/2 



r[a + f-f]r[6+|-f]r[a + 6+f-n]r[| 
T[a}T[b}T[e}T[a + b + e - n] 



(D4) 



3. d = 0, c / 0, e / and c = 0, d ^ 0, e/0 



For vanishing d with non- vanishing c and e we consider the integral H 2 {a,-, b, c, 0, e\n) 
for p 2 = m 2 , 

d n k l d n k 2 1 



H^\a,b,c,d,e\n) 



(27r) 2n [k 2 - m 2 + i0+] a [k 2 - m 2 + i0+} b [2p ■ k x + z0+] c 
1 



[2p ■ k 2 + i0+] d [2p ■k 1 + 2p-k 2 + i0+] e 



(D5) 



p 2 =m 2 N 



Note that the mass terms m in the first two propagators stem from the rescaling of the loop 
momenta, while we have to consider p 2 = m 2 N when evaluating the nucleon mass. In the 
calculations performed in this work the difference between p 2 = m 2 and p 2 = m 2 N in these 
integrals is of higher order. 

The result for H 2 1 ' 1 \a, 6, c, 0, e\n) is given by the sum 

c-l 



H^'^^a, b, c, 0, e\n) 



E 

1=0 



c-l 
I 



■l) l Z( c+l _ 2)/2 (a,b, c,0,e\n), 



(D6) 



where 



,-2-2a-2b-2c-2e 



Z a (a, b, c, 0, e\n 



m 



2n-2a-2b-2c-2e 



T[a + 1] 



x 



(4vr) n r[a + 2] 

r[a+j+_| - | ]Tp> + f + | - f]T[| + f]r[a + & + f + f - n] 
r[a]r[6]r[c]r[e]r[a + & + c + e-n] 
1, c/2 + e/2, a + c/2 + e/2 - n/2 
a + 2,a + 6 + c + e — n 



X 3 ^2 



(D7) 



d,e 



z) is a hypergeometric function. The case c = 0, d 7^ 0, e 7^ is obtained by 



and 3 F 2 

replacing c with d and interchanging a and 6 in Eq. (ID 71) . 
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4. c ^0, d^O, e^O 



For the case that none of the exponents c, d, e vanishes, it is convenient to perform an 
expansion into partial fractions, 



[2p ■ k 2 + i0+) [2p ■k 1 + 2p-k 2 + «0 H 
1 



[2p ■ ki + 20+] [2p ■ k 2 + i0+] [2p ■ ki + z0+] [2p ■k 1 + 2p-k 2 + i0+] ' 

(D8) 

until one obtains a sum of integrals of the form iff' 13 (a, b, c, 0, e) and H^' 1 \a, b, c, d, 0), 
which are evaluated as described above. 



5. Subtraction terms 

In addition to the integrals given above the evaluation of the subintegrals for the _F 4 terms 
requires the integrals 



iff' 13 (a, 6; win) = J 



d n k 



1 



(2vr) n [k 2 - m 2 + i0 + } a [2p ■ k + u + i0 + } b ' 



(D9) 



where u = 2p- q with g the second loop momentum. The integral if ] L ' 1 (a, 6; u\n) is given by 
the sum 



H[ 1 ' 1 \a, b; uj\n) 



l-2a-2b^n-2a-2b 00 r\ b _L 1 



m 
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(D10) 



The sum contains an infinite number of terms. However, when performing the second loop 
integration over q in the considered counterterm integrals, increasing orders of uo = 2p ■ q 
contribute to increasing chiral orders. Therefore only a finite number of terms in Eq. (IDlOj) 
is needed in the calculation of the nucleon mass. 



6. Results for iff' 13 (a, 6, c, d, e|n) and counterterm integrals 



The results for the iff' parts of the two-loop integrals contributing to the nucleon mass 
evaluated on-mass-shell are given by 
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where ft is the 't Hooft parameter and we have used ft 



1 E ~ ^ 

/_ — e^ - . In the above notation 



(An) 1 / 2 



the MS scheme of ChPT corresponds to subtracting all terms proportional to e 1 and setting 
fi = m. 

The corresponding counterterm integrals H^^(a,b,c,d,e) and H^^(a,b,c,d,e) are 
given by 
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The results for the products of one-loop integrals read 
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The integral H£ A \l, 0, 0, 0, l\n) can be written as H^'^il, 0, 0, 1, 0|n) by the substitution 
A; 2 i — ^ A; 2 + fci, and 

if^l, 0, 0, 1, 0|n) = ^(l, 0|n)^ (1) (0, l|n) = 0, 

since the infrared sing ular part iff^OlIn) of #i(01|n) = -U N vanishes. 

The counterterm integrals corresponding to the products of one-loop integrals read 
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